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Aerothermodynamics  of  Pulsed  Detonation  Engines 


Shmuel  Eidelman*,  Dmitri  Sharov*,  and  David  Book*^ 
*SAIC  McLean  VA  22102  /  '^nigmatics  Inc.,  Wash.  DC 

ABSTRACT 


With  no  moving  parts  in  the  power  production  section  and  a  thermodynamically 
efficient  cycle,  the  Pulsed  Detonation  Engine  (PDE)  offers  a  low-cost  alternative  to 
turbojet  and  liquid-propellant  rocket  engines.  Both  air-breathing  and  pure  rocket  modes 
of  engine  operation  offer  substantial  system,  materials,  and  cycle  advantages  that  will 
allow  us  to  bypass  the  scalabUity,  operational  range,  efficiency,  and  cost  limitations  of 
existing  engines.  We  introduced  the  modem  PDE  concept  in  1986^  and  experimentally 
demonstrated  its  operation  in  1986^  1994^  and  1998^  The  concept  was  subsequently 
treated  extensively  in  a  number  of  analytical  and  experimental  studies.^'*® 

This  paper  examines  the  thermal  balance  of  the  PDE  as  a  function  of  engine 
length,  cycle  frequency,  external  flow  Mach  number  and  other  parameters.  The  results 
indicate  that  PDEs  will  have  thermal  characteristics  that  are  drastically  different  from 
conventional  engines.  The  results  of  the  study  have  important  implication  for  PDE  design 
andanalysis. 

Introduction 

Thermal  management  is  a  subject  of  active  discussion  in  PDE  research.  One  point 
of  view  is  that  PDEs  will  allow  the  use  of  high-energy-density  fuels  without  a  need  for 
active  cooling  or  high-cost  materials  for  thermal  protection,  since  the  combustion 
chamber  walls  are  intermittently  exposed  to  hot  detonation  products  and  the  cold  fuel/air 
mixture,  with  the  cold  mixture  resident  longer  in  the  detonation  chamber  and  therefore 
playing  a  bigger  role  in  the  PDE  thermal  balance.  Another  point  of  view  is  that  the 
detonation  chamber  will  overheat  as  a  result  of  very  intense  heat  exchange  with  the 
imsteady  flow  of  detonation  products,  which  for  hydrogen/oxygen  detonation  can  reach 
temperatures  in  excess  of  3500  K.  It  is  known  that  some  groups  developing  PDEs  have  a 
problem  keeping  the  engine  walls  cool. 

Because  of  the  high  temperature  of  the  detonation  products  and  high  heat  transfer 
rates  from  the  nonsteady  flow  of  hot  gases  to  the  detonation  chamber  walls,  detonation 
chamber  cooling  is  a  major  design' consideration.  Analysis  of  the  heat  transfer  is 
complicated  by  the  fact  that  during  the  purging  and  refill  process  the  detonation  chamber 
will  be  fiUed  with  cold  gases  and  the  heat  flux  will  reverse  direction  from  the  hot  wall  to 
the  gas.  In  the  nonsteady  processes  of  fill,  detonation,  expansion  following  detonation, 
and  purge  of  detonation  products,  depending  on  engine  design,  different  parts  of  the 
engine  will  be  exposed  to  different  intermittent  flow  regimes.  The  problem  of  heat 
exchange  is  compounded  by  the  possible  problem  of  fuel  pre-ignition  that  can  occur  when 
the  fuel/oxidizer  mixture  is  injected  into  a  very  hot  detonation  chamber.  This  can  result  in 
combustion  of  the  injected  mixture  instead  of  detonation  and  in  some  cases  can  convert 
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the  chamber  into  a  flame  holder  that  will  make  subsequent  detonations  impossible, 
causing  complete  loss  of  thrust. 

hi  this  paper  we  consider  the  basic  dynamics  of  the  PDE  heat  exchange  for  a 
simple  engine  configuration. 


Problem  Description 


We  assume  that  the  PDE  device  has  the  simplified  configuration  shown  in  Figure 
1.  The  detonation  engine  consists  of  a  steel  tube  with  6-cm  diameter  and  2-mm  thickness 
(in  our  analysis  we  will  consider  various  tube  lengths).  The  air  is  injected  through  the 
forward-facing  inlets  and  mixes  with  fuel  injected  at  the  thrust  wall.  We  assume  that 
mixing  is  very  efficient  and  that  for  the  purposes  of  our  analysis  the  gaseous  fuel/air 


Purge  Phase 


Fuel  Injection  / 
Mixing  Phase 


Ignition  Phase 


Detonation  Phase 


Figure  1.  An  illustration  of  the  basic  phases  of  an 
airbreathing  PDE  cycle. 


mixture  fills  the  detonation  chamber  uniformly.  Detonation  is  initiated  at  the  thrust  wall 
by  a  thin  layer  of  high-pressure  gas  and  propagates  through  the  mixture.  This  is  only  a 
general  description  of  the  problem.  Below  in  our  description  of  the  numerical  simulations 
we  will  give  more  details  about  detonation  of  the  injected  mixture  and  the  detonation 
wave  parameters. 

The  fuel/air  mixture  has  the  temperature  of  the  ambient  air  and  fills  the  detonation 
chamber  gradually,  with  a  filling  rate  such  that  detonation  is  initiated  at  the  instant  when 
the  mixture  reaches  the  open  end  of  the  detonation  chamber.  The  cycle  period  consist  of 
the  time  needed  to  fill  the  detonation  chamber  at  a  certain  rate  and  the  time  it  takes  for  the 
detonation  wave  to  propagate  through  the  mixture  and  for  the  detonation  products  to  fully 
expand,  producing  thrust.  Thus,  for  example,  if  the  detonation  engine  is  1  m  long  and 
operates  at  100  Hz  and  we  know  that  the  detonation  and  expansion  processes  take  4  msec 
we  assume  that  the  fuel/air  mixture  is  injected  into  the  chamber  at  167  m/sec  to  allow 
complete  filling  of  the  detonation  chamber  in  6  msec. 
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Numerical  Simulations 

It  is  feasible  to  find  the  heat  flux  to  the  wall  during  the  process  of  injection, 
mixing,  detonation,  and  detonation  product  expansion  through  numerical  solution  of  the 
time-dependent  Navier-Stokes  equations  describing  viscous,  reactive,  and  turbulent  flow 
regimes.  Calculated  heat  fluxes  to  the  wall  can  be  used  as  boundary  conditions  for  the 
heat  transfer  problem  through  the  solid  walls  of  the  PDE.  Such  a  computation  can  be  very 
time  and  resource-consuming  and  at  the  end  its  results  cannot  be  fully  trasted  and  would 
need  experimental  validation  because  of  the  ambiguities  related  to  selection  of  the 
turbulence  model  and  spatial  resolution  requirements  needed  to  resolve  boundary  and 
mixing  layers.  Although  direct  Navier-Stokes  simulations  of  the  heat  flux  in  PDE  devices 
should  be  pursued  in  the  future,  in  this  study  we  adopted  a  simpler  approach  that  allows 
rapid  evaluation  of  the  PDE  wall  temperature  using  PC-based  simulations. 

In  our  study  we  use  solutions  of  the  Euler  equations  to  simulate  detonation  and 
expansion  processes  using  a  methodology  that  we  previously  developed  for  simulations 
of  these  processes  in  PDEs  (see  References  2-10  for  details).  We  assumed  that  the  PDE 
has  the  simple  cylindrical  geometry  shown  in  Figure  1,  where  the  tube  diameter  is  6  cm. 
We  also  assumed  that  the  inlets  are  closed  during  the  detonation  and  expansion  processes. 
The  detonable  mixture  of  stoichiometric  hydrogen/air  at  1  atm  and  300  K  is  initiated  by  a 
2  mm  vertical  layer  at  a  pressure  of  30  atm  located  at  the  thrust  wall.  The  numerical 
simulation  gives  us  the  Wstory  of  the  gas-dynamic  parameters  such  as  temperature, 
velocity,  and  density  on  both  internal  and  external  surfaces  of  the  engine. 

In  Figures  2,  3,  and  4  we  show  the  results  of  simulating  detonation  and  expansion 
processes  in  the  form  of  pressing,  velocity,  and  temperature  time  histories  for  a  17  cm- 
long  PDE  configuration.  We  assumed  that  the  engine  operates  at  100  Hz.  The  injection  of 
the  cold  mixture  begins  1  msec  after  initiation.  As  we  see  from  the  pressure  plots,  1  msec 
after  detonation  the  pressure  inside  the  chamber  has  the  ambient  value.  It  is  interesting  to 
observe  that  during  the  detonation  and  expansion  processes  all  points  along  the 
detonation  chamber  experience  the  same  pressure  and  velocity  conditions  (except  during 
a  short  period  when  the  detonation  wave  propagates  and  transition  to  the  C-J  condition 
takes  place).  However,  the  temperature  regimes  experienced  by  the  locations  along  the 
wall  are  a  very  strong  function  of  location  as  is  evident  from  Figure  4.  The  point  2  cm 
from  the  thrust  wall  will  be  at  a  high  temperature  of  -2000  K  only  2  msec  and  at  300  K 
for  8  msec.  At  the  same  time,  the  point  at  10  cm  from  the  wall  will  be  exposed  7  msec  to 
~2000  K  and  3  msec  to  300  K.  The  points  near  the  end  of  the  detonation  tube  will  have  a 
lower  temperature  at  the  end  of  the  exp^sion  cycle  due  to  suction  of  the  ambient  air  after 
the  initial  overexpansion  of  the  detonation  products. 

We  assume  that  the  PDE  device  is  flying  with  a  given  flight  Mach  number  at  sea- 
level  conditions,  where  the  external  flow  is  cooling  the  detonation  chamber  wall.  The 
parameters  obtained  as  a  result  of  simulation  of  the  detonation  and  expansion  processes 
and  the  parameters  of  the  external  flow  over  the  PDE  are  used  as  boundary  conditions  for 
the  heat  transfer  problem,  which  is  solved  for  every  point  along  the  engine  waU. 
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100Hz 


Time  (sec) 

Figure  2.  Pressure  time  history  at  different  locations  of  the  PDE  detonation 
chamber  as  function  of  cycle  time  for  a  17-cm-long  engine  operating  at  100  Hz. 


100Hz 


Time  (sec) 


Figure  3.  Velocity  time  history  at  different  locations  of  the  PDE  detonation 
chamber  as  function  of  cycle  time  for  a  17*cm*long  engine  operating  at  100  Hz. 


4 


PAGE 


JAN-12-00 


15:14  FROM:SAlC/ATG  McLean  VA 

Presented  at  36*  lANNAl  ;  ,?.  Meeting,  Cocoa  Beach,  FL,  1999 


ID:+1  703  G76  5323 

'  w  - 


7/1 


Figure  4.  Temperature  time  history  at  different  location  of  the  PDE  detonation 
chamber  as  function  of  cycle  time  for  a  17-cm-long  engine  operating  at  100  Hz. 


The  heat  transfer  problem  is  described  by  the  heat  transfer  equation  in  cylindrical 
coordinates: 


ar 

dt 


=  K 


ar^  ^  r  ar 


The  boundary  condition  at  r=Ri  is 


^  (t  _ 7’ 'j 


The  boundary  condition  at  r=l?2  is 


h  = 


-r ) 
ar 

Nu--k 


K  =  - 


where  h  is  the  heat  transfer  coefficient  for  the  gas;  fcw  is  the  thermal  conductivity  for  the 
wall;  k  is  the  thermal  conductivity  for  the  gas;  and  I  is  the  length  scale.  In  our 

calculations  we  used  JC=12.5xlO'^m^/s  (steel),  fcw  =45.4  W/(m  K)  and 

k  =  Pr ,  where  Pr  =0.72  and  y=1.4 
y-1 

The  viscosity  is  computed  from  the  Sutherland  law. 
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The  Nusselt  number  on  the  external  surface  of  the  engine  is  calculated  using  the 
formula  for  turbulent  flow  over  a  flat  plate,  which  according  to  Reference  1 1  is  also 
applicable  for  external  flow  over  a  cylinder; 

iVn=0.03Re"-®Pr°"^ 

The  scaled  length  used  for  calculating  the  Reynolds  number  is  equal  to  the  distance  from 
the  thrust  wall. 

The  Nusselt  number  on  the  internal  surface  of  the  engine  is  calculated  using  the 
turbulent  flow  in  the  tube  approximation: 

/  iVM  =  0.021Re®^Pr®''' 

Here  the  length  scale  used  for  the  Reynolds  number  calculation  is  equal  to  the  internal 
diameter  of  the  tube. 

In  all  calculations  we  assumed  that  the  next  injection  is  started  1ms  after  the 
ignition  of  the  previous  portion  of  the  detonable  mixture.  As  explained  above,  the  speed 
of  the  fuel/air  injections  for  different  cases  was  adjusted  to  allow  complete  filling  of  the 
detonation  chamber  volume  for  a  given  cycle  period.  In  Table  1  we  list  the  injection 
speeds  that  were  used  in  different  simulations. 


Table  I.  Injection  speed  (m/sec)  for  different  engine  length  and  operation  frequency. 


Engine 

Engine  operation  frequency 

length 

SOEz  \ 

lOOEz 

lOOBz 

400Hz 

17cm 

7.5  m/s 

16.7  m/s 

37.5  m/s 

100.0  m/s 

34cm 

15.8  m/s 

33.3  m/s 

75.0  m/s 

— 

64cm 

31.5  m/s 

66.6  m/s 

150.0  m/s 

— 

Results 

Our  methodology  allowed  us  to  simulate  the  PDE  wall  temperature  as  a  function 
of  time,  engine  length,  external  flow  Mach  number,  and  engine  operation  frequency.  To 
examine  the  effect  of  only  the  four  parameters  listed  above  we  performed  aU  the 
simulations  for  a  hydrogen/air  detonation  in  a  6-cm-diameter  engine;  however,  this  is  not 
a  limitation  of  OUT  methodology,  which  can  be  used  to  calculate  thermal  balance  for 
complex  geometries  and  other  detonable  mixtures. 

In  Figure  5  we  show  engine  wall  temperature  as  a  function  of  time  and  distance 
from  the  thrust  wall  for  a  34-cm-long  engine  operating  at  100  Hz  in  an  M=0.3  external 
flow.  We  see  in  this  picture  that  the  wall  temperature  will  reach  steady-state  values  only 
afteF.6Qsec  of  engine  operation  (or  6000  cycles).  Note  that  we  assumed  a  wall  thickness 
^^(pf  2nm^  for  thicker  walls  or  lower  heat  conductivity  the  time  required  to  stabilize  the 
tSnp^ture  can  be  even  longer.  This  property  implies  that  the  heat  losses  of  the  PDE 
cycle  are  low,  which  is  another  indication  of  cycle  efficiency.  The  relatively  slow  heating 
.  of  the  wall  can  be  exploited  in  engine  design,  particularly  for  engines  with  short  operating 
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times.  By  choosing  materials  with  a  low  heat  transfer  coefficient  one  can  operate  in 
transient  mode  for  minutes. 


100Hz,  34cm  engine 


Figure  5.  Engine  wall  temperature  as  a  function  of  time  and  distance  from  the 
thrust  wall  for  a  34  cm  long  engine  operating  at  100  Hz  in  M=0.3  external  flow. 


100Hz,  M=0.05, 34cm  engine 


Figure  6.  Engine  wall  temperature  as  a  hmction  of  time  and  distance  from 
the  thrust  wall  for  a  34  cm  long  engine  operating  at  100  Hz  in  M=0.05 
external  flow. 


7 


PAGE  10/1 


-  JAH-.12-00 


15:16  FROM:SAIC/ATG  McLean  VA 

Presented  at  36*  JANNA^^ceting,  Cocoa  Beach,  FL,  1999 


ID:+1  703  676  5323 

c 


Another  property  ,  thk  is  clearly  identifiable  .from  Figure  5  is  that  the  wall 
temperature  is  a  strong  function  of  the  axial  distance  from  the  thrust  wall.  -  The  engine 
wall  temperature  2  cm  from  the  thrust  wall  is  320  K,  while  that  at  the  aft  end  the  wall 
reaches  710  K.  This  can  be  explained  by  the  different  time  of  exposure  to  hot  and  cold 
gas  at  these  locations.  The  section  of  the  engine  2  cm  from  the  wall  be  exposed  most  of 
the  time  to  the  cold  flow  of  injected  mixture,  while  the  point  at  the  aft  end  will  be 
exposed  to  the  cold  mixture  for  a  very  short  time. 

In  Figure  6  we  show_wall  temperatures  for. the  same  engine. .as,  in  Figure  5, 
operating  in  a  Mach  =  0.05  external  flow.  Here  we  see  that  the  engine  waU  temperature 
can  reach  1200  K.  However,  it  is  interesting  to  note  that  the  wall  heats  up  significantly 
more  slowly  than  in  previous  case.  It  takes  twice  as  long  to  reach  maximum  temperature 
in  the  M=0.05  external  flow  as  in' the  Af  =  0.3  external  flow. 


engine  operating  at  Mach  =  0  J  external  flow  conditions. 


In  Figure  7  we  examine  the  maximum  wall  temperature  at  different  locations  as  a 
function  of  engine  operating  frequency  for  a  17-cm-long  engine  operating  in  M=0.3 
external  flow  conditions.  This  figure  Dlustrafes  that  the  wall  temperature  is  a  linear 
function  of  frequency..  The  wall  temperature  increase  is  almost  exactly  proportional  to 
frequency  increase.  Thus,  an  eight-fold  increase  in  frequency  (and  power  density)  from  50 
Hz  to  400  Hz  leads  to  about  an  eight-fold  increase  in  the  wall  temperature  difference  and 
ambient  temperature.  This  is  frue  for  all  locations,  but  the  wall  section  near  the  aft  end 
where  the  effects  of  overexpansion  and  suction  of  the  ambient  air  make  this  relation 
between  the  temperature  and  frequency  more  complicated. 

.  In  Figure.  8  we  present  a  summary  of  our.  simulation  results  for  17-cm,  34-cm,  and 
64-cm-long  PDFs  operating  at  50  to  200  Hz.  All  simulations  where  conducted  at  Mach  = 
0.3  flight  conditions.  This  figure  shows  that  longer  engines  will  heat  up  substantially  even 
when  operated  at  a  relatively  low  100  Hz  frequency.  The  maximum  temperature  reached 
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in  the  aft  section  is  a  ftinction  of  the  thrust  generated  by  the  engine.  For  example,  a  34-cm 
engine  operating  at  200  Hz  will  produce  about  the  same  thrust  as  a  64-cm  engine 
operating  at  100  Hz,  and  they  will  have  about  same  maximmn  temperature  at  the  aft 
section.  The  same  is  true  for  a  17-cm  device  operatiag  at  200  Hz  and  a  64-cm  device 
operating  at  50  Hz. 


Figure  8  Maxunum  wall  temperatures  as  a  function  of  length  for  PDFs 
operating  at  different  firequencies  and  at  Mach=0.3  external  flow  conditions. 


Conclusions 

We  have  developed  a  methodology  that  allows  numerical  simulation  of  the  PDE 
wall  temperature  and  heat  losses.  Using  this  methodology  we  examined  the  main  trends 
of  detonation  chamber  wail  temperature  dependence  on  engine  length,  frequency,,  external 
Mach  number,  and  time.  Although  demonstrated  for  simple  cylindrical  configuration,  this 
methodology  is  general  and  it  can  be  used  for  evaluation  of  the  waU  temperatures  of 
complex  3-D  PDE  designs.  As  a  result  of  our  study  we  can  draw  the  following 
conclusions: 

•  Due  to  the  intermittent  exposure  to  cold  and  hot  gas  the  PDE  wall  will  heat  up 
relatively  slowly  reaching  maximmn  values  for  the  simulated  configuration  in  1  to  2 
minutes  after  engine  start. 
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•  The  PDE  will  reach  its  final  temperature  more  quickly  for  higher  subsonic  flight 
Mach  numbers; 

•  For  the  chosen  cylindrical  configuration,  the  wall  temjjerature  will  increase  as  a 
linear  function  of  distance  from  the  thrust  waU,  and  the  sections  closer  to  aft  end  will 
be  the  hottest; 

•  Longer  engines  will  reach  higher  maximum  wall  temperatures  for  the  same  engine 
operating  frequency; 

•  Wall  temperature  increases  linearly  as  the  engine  frequency  increases; 

•  For  the  chosen  PDE  configurations  the  maxiinnTn  engine  temperature  is  a  fimetion  of 
engine  thrust  regardless  of  the  frequency  and  length  combination  that  allows  us  to 
obtain  the  thrust  value. 

The  new  simulation  methodology  can  be  used  to  design'PDE  configurations  that 

win  not  need  forced  cooling,  thus  significantly  simplifying  engine  design. 
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